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relays AMIS 28 


This is a report of a seareh for propagation anomalies using a 
large quantity of high frequeney data produeed as a byproduct of 
BRIGHAM, a Department of Defense projeet. The BRIGHAM data is 
based on 890 KHz wide samples of the HF’ spectrum at a 25 eycle rate, 
Wsingee 2,6 KHz resolution for a duration of approximately 2.5 minutes. 
This method of data eolleetion is unique and it was hoped that propaga- 
tion anomalies, ineluding wide band anomalies, might be detected. 
Anomalies are believed to occur in the propagation of radio signais 
and they are usually other than known, routine effeets but may include 
known effeets which cannot adequately be explained. The seope of this 
examination was limited to the visual analysis of eomputer proeessed 


data presented on an interaetive graphies unit. 
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i OBGEGLIVES 


his’ project provided thie opportunity to cxabmine unique data 
gathered with the AN/FLR-15 wide-band scanning receiver and the 
AN/FRD-10 circularly disposed antenna systém. The project was 
suggested by the availability of this large volume of data which had 
been gathered at great expense for another purpose. The data seemed 
fomoltemeereat potential andthe possibile basys for considierabletnesicarch 
and understanding about electromagnetic propagation. The initial con- 
cept was to observe data in the different form of a 3-dimensional con- 
uMmovani plitude wirequichnay and time, This concept included the use 
of a computer and interactive graphics unit to vary the size and dim- 
ensions of the contour to obtain a new view of the high frequency region 
of the spectrum and perhaps observe propagation characteristics not 
easily recognized otherwise. 

The objectives of this project were to take the given sets of 
data, write the computer programs necessary to present the data 
visually and to analyze the data. The objectives of the analysis were 
twofold: To study propagation effects on high frequency signals, 
becoming familar with the omy the mode of visual presentation 
and the appearance of routine propagation effects; To look for unusual 


and anomalous propagation effects. 





foeeN SRO C LION 


It is stressed that the scope of this effort was confined toa 
visual analysis of the data for detecting unusual patterns and anomalies 
which could be attributed to propagation conditions. The search 
was conducted in the high frequency region (3-30 MHz) using data 
on computer cards. The data cards were produced during a phase 
Gi RiOHA Moe@amproject conducted for the Department of Defense by 
Sanders Associates, a private corporation. Each data set represents 
the output of a swept-tuned panoramic receiver and preserves amp- 
litude information vs frequency and time. For visual presentation, 
the data were processed ona Scientific Data Systems XDS 9306 
computer and displayed onan Adage AGT/10 graphics unit (figure 1). 
The equipment is part of the Electrical Engineering Computer Center 
of the Naval Postgraduate School. Interesting sets of sweeps have 
been extracted as one picture, plotted on a CALCOMP 563 plotter 
and included in this report. Forty-six data sets were available for 
examination, cach consisting of 7195 cards. Several data sets were 
examined in detail with the intent of using other sets for verification 
of observed anomalies. Recommendations for further examination 


are given at the end of this report. 
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A, ee Ol rh DATA 

The data were gathered using a swept-tuned receiver and either 
a Omnidirectional or a 12° beam antenna. Each data set was produced 
by sampling the receiver output as it scanned downward through an 
890 KHz band. Sets were gathered at different base (starting) fre- 
Guencies amd al different times! of the day. Thc! slets! contain the 
amplitude (in db) of the 318 contiguous frequency bins (sampling points) 
and consist of 7194 cards to represent the 3597 scans. The initial 
card in each set contains date, time and parameter information (fre- 


Gacnieyer ete.) tor that set. 


Be Pio we iON OF THe DATA 

Computer programs for (1) transferring the data from cards to 
-tape, (2) visual presentation of the data and (3) extraction of the data 
for graph plotting are included at the end of this report. There are 
eve Computer proprams for presenting the data and they reflect the 
trade-off between the width of the displayed spectrum and the number 
of scans which could be viewed Simultaneously. This restriction was 
awe to core memory size in the XDS 9300 and again in the AG / VOR 

We Anomaly A 

This program was used for viewing the data as a contour 

and required appreciable dimensions in time and frequency to be of 


use. It was initially adapted froma library program but, along with 





two accompanying subroutines, was changed almost entirely in adapt- 
ing it to the BRIGHAM data. Anomaly Plot A was used with this 
program for plotting graphs. The data were displayed as shown in 
figure 2 with the most recent scan at the bottom of the picture. For 
each picture change, a new scan was brought on at the bottom and the 
oldest scan moved off at the top. The scans were not moved up in 

a continuous or movie-like fashion but all shifted upward a stepata 
time to replace previous scans. Options were available for bringing 
On more than one scanata time. The amount could be varied so as 
to change the entire picture at once if desired. Scans could be skipped, 
taking every fifth or tenth scan for example, effectively decreasing 
the receiver sampling rate to accentuate long term or slowly changing 
trends in the data. Also, the time axis could be rotated to vary the 
3-dimension or surface contour effect. Other features are described 
in the computer program section. An optimum size of 130 bins and 
20 scans was chosen for this program. A deficiency in this presenta- 
tion was that nde changes with time were not easily observed. 
eo t. Andee this, the Anomaly B program was written. The scale 
at the bottom of figure 3 and all the graphs is 10 frequency bins 

(28 KHz) per division. The leftmost bin in each group is labeled with 
the computer variable ISTRT, for the integer value of the ''starting" 
bin. ISTRT is a NAMELIST variable and its value can be changed 


at the graphics console. It enabled the viewer to look at different 
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Portions @t ime Bamd-) For identification, bins'#re™located by counting 
to the right from the starting bin (frequency decreasing toward the 
right). The vertical scale of the Anomaly A presentation is 5 db of 
amplitude per division (referred to receiver sensitivity). Scans are 
spaced 1/25 second apart in time or multiples thereof when scans are 
skipped. 
Ge Anomaly B 

This program was used to examine more cerefully the 
changes in amplitude among a number of signals (see figure 3). At 
the top are five traces of signals chosen by specifying the bin number. 
The traces show the amplhitude history of each signal up to a maximum 
strength of 30 db each, beginning with the present base scan at the 
left and extending to the right, showing the last 80 scans. The am- 
plitude was displayed after taking the maximum value of the signal 
in three bins (the one selected and one on either side) and scaling 
it. At the bottom are displayed all 318 bins to aid in selecting the 
five signals to be eee Because of the memory size trade-off, 
only five Sicans are ishown simultaneously. ‘Zhe vertical scales 
differ from those in Anomaly A. Because of a smaller scale, the 
divisions are 10 dbapart. Anomaly Plot B was used to plot this 


display on paper. 
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aX PRELIMINARY OBSERVATIONS 

Before beginning the visual analysis, it was useful to predict 
the type of displays the data could yield. These preliminary observa- 
tions helped visualize numerous features and signal variations not 
due to propagation so that they could be recognized and not be inter- 
preted as anomalies to be checked and verified. Examples of these 
features included the following: 

Ie leo We tomete mecha Gderet IStlCS. «5 @urning on and oft; dilierent 
types of keying which affect frequency and/or amplitude. 

a Receiver characteristics: Filters for notching out inter- 
fering signals. 

Other preliminary#observations! predicted limitations on the 
ability to observe valid propagation characteristics which might be 
present. These limitations were due to receiver characteristics, 
the nature of the data and the method of visual presentation. 

ie limitations of the receiver: The 890 KHz band restricted 
the ability to relate the effects on harmonics in the megahertz region 
and the ability to observe variations over a wide band. 

ae Pi totanvoiceauesto tne data mune lack of resolution 
imposed by the 2. 8 KHz wide frequency bins was expected to make 
Pe inapocsibic tomecasure slipnte variations imivequency. Better 
resolution could provide better identification of the types of signals 


present for the purpose of selecting ones with more stable transmitted 





characteristics to permit more reliable analysis and measurement. 
The short duration (2.5 minutes) sample provided a mountain of data 
but was expected to prevent observation of some propagation charac- 
teristics which require a longer period to manifest themselves. 

oe Shortcomings in the visual mode of presentation were not 
Poresieen until much later in the project. At first only a general con- 
Ceperot the type of presentation existed and the final characteristics 
were considerably different. In fact, Anomaly B was developed 
only after much analysis and frustration with Anomaly A. It should 
also be noted that although the B version was developed last, it was 


not a replacement and both programs complement each other. 


Is). PROCEDURE 

After reading a data set onto tape, the data were analyzed 
Visine the “nomaly A presentation to examine the entire length and 
width of the set. The object at this stage was to find gross anomalies 
iiPttewomdinaemshonal contour which mipht*te easily recognized. 
Hypothetical examples hoped for might have been any of the following: 
Peatiaseianve one ueelon only, of the comtour;isimilar changes to 
one or more amplitude humps (signals or combinations of signals) in 
separated but perhaps related areas of the contour; holes or the 
absience of signals in limited bands of the contour. The Anomaly B 


program was then used for more critical analysis of amplitudes over 


IZ 





@ longer period. ‘The Avgayslitendle xis save traces displayed data over a 
period four times as long as in Anomaly A (80 data points vs the 20 
scans) and examined it from a different perspective (abscissa in 
dimensions of time rather than frequency to examine individual 


signals). 


Cc. FINDINGS AND OBSERVATIONS 

Datatisetsi4, 5,°5, 9, 15, 32, 34 and 147 were examined in 
detail. Set 4 was the first set examined and received considerable 
attention while being usied to develop the first computer program. A 
routine for numbering each scan and later displaying the number of 
mac base Or bettom scan on the Screen was developed and became 
indispensable in helping identify the location of anomalies or unusual 
data for subsequent observations. 

The first unusual feature noted is shown in figure 4 (upper left 
hand corner). The amplitude hump there is about 4-6 KHz wide and 
useseen to move-downward in frequency with time at a rate of about 
250 KHz per second. An effort to trace it through other signals was 
unsuccessful. The left hand margin represents the upper limit of 
the band (J8 MHz), therefore tracing if at higher frequencies was not 
possible. The rate of frequency change was too fast to be classified 
as a ‘whistler’! or "dawn chorus''. These phenomena are only found 


at considerably lower frequencies (several kilohertz) and the bandwidth 
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FIGURE 4. Data set 4 


Base scan 2109 
Base frequency 18 Mliz 





of the amplitude hump was too great for this type of classification. 
This was the only occurrence of this phenomenon which was observed. 
It may have been a transmitter being tuned or a special purpose signal. 
Pelco ii.c pec entimiascis 4, 5, G99 15, 34 and 147 tare sudden 
disicontimuities across the entire spectrum of the set. Examplestare 
shown in figures 5 through 7. The transitions were too rapid to be 
attributed to propagation effects and a satisfactory explanation has 
not been found. Late in the analysis, in set 34, several pictures from 
Anomaly A were pieced together (figure 7). The figure shows an 
effective shift of 160 scans. At one point it seemed to be the result 
ofa ae data card. Two cards are required to describe one scan 
and if one is missing, the next card and all subsequent cards are all 
shifted back 160 scans. It would give the appearance on the left half 
of the scan being exchanged with the right half, as it appeayrs in 
figure 7. A second missing card would then put them back in the 
right order, ending the "discontinuity". Unfortunately, this simple 
explanation was not correct because an examination of the data in the 
region of the set 34 discontinuity revealed no missing or out-of-place 
cards. Examination of cards in other sets with discontinuities also 
revealed no missing cards. The bad data therefore must have been 
CUestOousOimie SOtipolnecciver or recording malfunction. 
Figure 8, instead of showing any anomalies, shows a set of 


well ordered, continuous signals of nearly constant amplitude and 
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FIGURE 5. 
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FIGURE 6G. 
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Sia vo meeioesmeonnrcastead ihn figure 9 by an irrepular contour. In 
addition to variations in amplitude, some signals seern to be turned 
on and off alternately. The region between bins 30 and 38 contains 
omsHemal rcoularly occurring only once in nearly every three scans. 
At bin 60, the pattern is irregular and is probably an on-off keyed 
signal. 

When figure 10 was first observed, it was thought that a genuine 
anomaly has finally been found. Band limited noise seemed to appear 
first to the right of bin 160 and then spread to higher frequencies 
before subsiding wlhics presience of the notchisdemed connected also, 
marking the lower limit during the first scan and persisting before, 
PonouUpwomr and attemtnie NO1IScun tiine. After some study and con-— 
sideration it was realized that each sweep required a finite amount 
of time and that the onset of noise appeared, although abruptly, at 
midscan. It is seen to grow weaker with time bcfore finally 
disappearing. 

marine: Prelysis was conducted with Anomaly B to examine 
amplitude variations of individual amplitude humps. Examination 
was begun in set 32 because it was free of the discontinuitics present 
in cach of the previous scts. Also, set 32 was takcn froma 12 beam 
antenna rather than the omnidirectional one used in the sets from 
4 +0 15. The narrow beam was chosen to confine the signals observed 


FO ad patticular repion of oriwin. Although points of signal origin, 
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and consequently propagation conditions, would still vary with each 
signal, the variations would not be as extreme as in the 360° case. 

Picwm@cs melewle and) lo were extracted from set 32 and compare 
the amplitudes of signals found in bins 19, 64. 101) 156 and 182. 
Each of thelthree figures were drawn with a different sampling rate 
but have the same terminal or base scan (leftmost points at the top 
of the figures and the bottom scan in the lower section of the figures). 
Figure 11 shows every other scan ‘lag variable NSKIP set equal to 1 
in the computer program) beginning with scan 2840, covering a 6.4 
second period. Figure 12 shows every tenth scan (NSKIP=10) begin- 
ning with scan 2200 (covers 32 seconds). Figure 13 shows every 35th 
scan beginning with scan 200 (covers 112 seconds). The scans at the 
bottom of each picture represent the entire data spectrum and show 
tac last five scans at the particular skip rate noted. 

A few similarities exist between the traces in figures 11 through 
13. Propagation effects had been expected to create corresponding 
dips in amplitude among the signals of similar points of origin (such 
as — SOW SOULILNew welte ExpeCuadio appCar over longer periods. 
On the other hand however, the rapidity of the fluctuations is not 
unusual. lJonospheric surface irregularities and layer displacements 
could account for them as described in reference 1. David and Voge 
[Ref. 2], for example, report that horizontal Meola Sees have been 


measured by doppler techniques at velocities up to 60-80 miles per 
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hour. At the top right in figure 11, traces 2 and 3 show an inverse 
relationship. Traces in figure 12 indicate that the higher frequencies 
were affected first with the effects shifting lower with time. It is 
realized that since the signals selected for observation here were 

not equally spaced in frequency, the straight line analysis may not 
be particularly valid. Also, because of the small number of similar- 
ities in these samples, they may only be coincidental. 

Figures 14 through 17 were also taken from set 32 but compare 
the amplitudes of bins 35, 85, 167, 224 and 303. The sampling rates 
are shown for each figure as wellas the terminal scan number and 
duration of the sample. Figure 14 shows no anomalies or correlation 
among signals with NSKIP=2. In figure 15, a periodic relation of longer 
duration seems to exist between traces 1 and 3. The propagation con- 
ditions affecting the two signals are similar but with the higher fre- 
quency signal (tuvccelaat the top) being affected earlier in time. These 
effects are seen again in figure 16 but at a higher skip rate (dips 
closer together) Oyere they can be observed over a longer period. 

No At exists to verify this pattern for the preceding time interval 
butit does not appear in the following interval (figure 17 - covers 

the interval from scan 1600 to number 3200 ~ NSKIP=20 over 80 data 
points). It may have died out or the display characteristics may have 
prevented it from being seen. The relationship stands out because 

of the frequency dependence of the effects. This type of anomaly is 


characteristic of what was expected when the project was begun. 
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Figure 18 shows an unusual relationship between the signals of 
traces 1, 2 and 3. Two of the signals experienced some sort of. 
attenuation while the signal in the middle was enhanced. A similar 
inverse relationship was observed in figure 1]. The same phenomenon 
Poemmlcugesho Can be seen in figure I but at a higher skip rate. Also 
present were periodic variations in amplitude, particularly in trace 3. 
Periodic variations became more pronounced in trace 4 (figure 20) 
in a later time interval. 

In addition to observing those groups of signals already noted, 
other groups of signals were also chosen in order to sample a wide 
range of possible combinations. Some groups spread the sample over 
the entire data spectrum while others were confined to a narrow range 
of bins. A spread of low and high skip numbers was also included. 
Specifically, the following groups of signals in set 32 were among 
those examined but with nothing significant observed: Bins 3, 9, 16, 
2deand 2omaith NoKIP equalto 1,747 16 and 40; Bins'135, 147, 158, 

168 and 177 with NSKIP Cello tem mele aides bins Oye bs, 17,922 and 


27 with NSKIP equal to 2, 10 and 35. 
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ive CONCLUSIONS AND RECOMMENDATIONS 


Although a tremendous amount of data was available as a by- 
product of project BRIGHAM and a large quantity of it was examined in 
detail during this effort, the many anomalies hoped for at the begin- 
ning were not seen. Three significant phenomena were observed and 
Poesetved in this report. The first, the signa! decreasing in freq- 
uency with time was almost certainly not a product of a propagation 
anomaly. Its identification might prove to be an interesting challenge, 
however. The second phenomenon showed that its periodic features 
were frequency dependent. Although periodic fading or intensity 
fluctuations were noticed at various times while examining the data, 
none were so well defined nor were they so well related between 
signals as inthis case. Comparing the ionospheric layers to a cloud 
whose features are constantly changing, it was not surprising to see 
Perc prcnomenonm contined to a Short period of time, Likewise, the 
irregular features of the ionosphere would allow several patterns 
to occur at once, superimposed on each other to create the overall 
observed effect. Some of these secondary features seem to show 
through in figure 16. The third phenomenon was the simultaneous 
enhancement and attenuation of signals. The seemingly incompatible 


conditions which caused this may have only been due to multipath 
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effects. It is strange, however, that no such pronounced effects 
ea observed before or after the anomaly, particularly in trace 3, 
and that they occurred simultaneously among the three signals. It is 
expected that a great number of other anomalies exist and that many 
of them are preserved in the BRIGHAM data. Approximately one 
fifth of the available data were examined and due to time limitations, 
only one fifth of that amount was examined with Anomaly B. Analysis 
with Anomaly B was particularly time consuming because only five 
signals could be examined at once. It is likely that other anomalies 
are observable with the existing programs and many others detectable 
with modifications and different analysis tools. Some recommendations 
for improving the analysis techniques are as follows: 

ts Amend the Tape Load program to read the binary number 
of each data card and alert the operator to missing cards. A one-up 
numbering system to keep pace with card numbers could be used 
with an ''if"' statement to trigger an output statement when a card is 
discovered missing. | 

ce Modify the procedure for selecting data for the amplitude 
history traces in Anomaly B. When cards are skipped to examine 
long-term propagation cffects, an averaging routine should be entered 
to remove the short term amplitude fluctuations. The presence of 
the short term cffects tends to mask those manifested only over a 


longer period of time. A modification of this type was attempted 





near the end of the analysis but was unsuccessful due toa lack of 
time needed to perfect it. 

or Incorporate a feature for backing up the tape by a variable 
(NAMELIST variable) number of scans. Two attempts to do this 
failed. In one case it was due to the fact that each scan constituted 
one ''physical record'' in length and that backspacing was based on 
"logical records'', If accomplished, the feature would be a significant 
time saving addition. 

4. Temporarily abandon the visual analysis approach of 
Anomaly B and employ statistical analysis techniques. Evaluation 
of the correlation between various sets of signals could be conducted 
rapidly and rather exhaustively on a larger computer such as the 
IBM 360 available at the Naval Postgraduate SeWoqime bare tli nan 
analysis with Anomaly A might prove useful in the selection of signals 
for analysis and indispensable in locating the discontinuities described 
earlier so as to avoid misinterpreting results. 

Ds Use of a larger, faster computer suchas the IBM 360 to 
load a onto magnetic tape would significantly reduce "dead time" 


in an extension of this project. Over an hour was required to loada 


single data set onto tape with the equipment described in this report. 





GIO Sere Ne 


i. Bin (frequency bin): A 2.8 KHz wide data sampling point. 
Energy present in these limits was quantized and represented in the 
pan eeptrom 0-25 db. 


a Data set: 3597 successive sampled sweeps or scans of the 
receiver and preserved on punched cards of magnetic tape. 


Se Dawn chorus or whistler: Interference in the audio range from 
storm discharges or high energy particle interaction in the ionosphere. 
Dawn chorus is produced around sunrise and sounds like a serenade 
of birds. lft is a rise in frequency above about 2 KHz. Whistlers are 
audible decreases in frequency and are propagated from storm dis- 
charges along magnetic force lines of the earth in the range froma 
few hundred hertz to 20-30 KHz. 


4. NAMELIST variables: Computer variables, properly specified, 
which can be entered or changed after the program has been compiled 
and is running. These variables were used to control various aspects 
of the program. The most used were: 


1S eae Inverter NUmMbe ceomthe Veitmost or Starting binjto be 
displayed. 

Petite cer limber olscddiacem bins to be displayed. 

INS SG Nev ee Number of scans to be skipped over. 


NMBRJIMP ~ Number of scans to be brought on (and moved off) 
the screen at one time. 


ae Bical neeswee Qo: the mecerver irom one end of the band to the 
uther. 
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The following subroutines, GINP, GINPUT, FNS and VCD, are 
common to and required by both programs Anomaly A and Anomaly B. 


They are library subroutines but are nonstandard and subject to change. 
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"GRAPH PLOT" is®a library program on paper tape and 1s\ usted with 
the CDC 160 computer to drive the CALCOMP 563 plotter. It is 


found in the Electrical Engineering Computer Center of the U. 5S. 


Naval Postgraduate School. 
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